Boll. Malacol., 43 (9-12): 131-138 (2007) 


Canonical morphometry versus statistical treatment of outlines 
through Fourier shape analysis: an empirical comparison 


Francesco Brusoni* & Daniela Basso* (El) 


* Via S. Maria del Sole 10, 
26900 Lodi, Italy, 
f.brusoni@katamail.com 

# Dip.to di Sc. 

Geologiche 

e Geotecnologie, 
University degli Studi 
di Milano-Bicocca, 

Piazza della Scienza 4, 
20126 Milano, Italy, 
daniela.basso@unimib.it, 
(s) Corresponding 
Author 


Abstract 

The analysis of outlines is considered a powerful tool for morphometric studies, since the absence of assump¬ 
tions on characters to be selected and measured reduces subjectivity. Fourier shape analysis converts each di¬ 
gitalized outline into a suite of numerical descriptors (the Fourier coefficients) that describe curves. The coeffi¬ 
cients are the input to multivariate statistical analyses, that elaborates a numerical comparison of shapes. A 
friendly, easily accessible and practical method for a statistical approach to shape analysis and the sensitivity 
of both Elliptic Fourier Analysis and Fast Fourier Transform have been tested here on a set of Goodallia trian¬ 
gularis shells. Shape analysis results to be faster and more reliable than canonical morphometry. In the com¬ 
parison of the Principal Component Analyses based on the two different Fourier transformations, Fast Fourier 
Transform produces the sharpest separation of the shell outlines into clusters while Elliptic Fourier Analysis re¬ 
sults are comparably less clear. Both techniques easily separate the two forms of Goodallia triangularis. 

Riassunto 

La descrizione quantitativa delle forme biologiche viene tradizionalmente affidata alia morfometria. Nella 
routine, la tecnica morfometrica viene tuttavia applicata a un limitato numero di caratteri misurabili, scelti 
in modo piu o meno arbitrario, riconducendo quindi la variability registrata a pochi parametri, Molte carat- 
teristiche morfologiche vengono tradizionalmente descritte in modo qualitativo (ad esempio la convessita 
della spira nei gasteropodi), senza quindi assegnare alcun parametro matematico di misura che possa suc- 
cessivamente consentire un confronto su base statistica-quantitativa. 

Per evitare che la scelta arbitraria di alcuni caratteri rispetto ad altri influenzi il risultato dello studio morfo- 
logico-comparativo sarebbe necessario considerare la totality dei caratteri misurabili. II tentativo di avvici- 
narsi a tale situazione ideale, evitando una eccessiva semplificazione, richiede tuttavia un enorme dispen- 
dio di tempo. A questo si aggiunga la percentuale di errore dovuta alia misurazione manuale e I'aumento 
delle difficolta quando le dimensioni degli oggetti di studio sono molto ridotte. 

La trasformazione del profilo di una struttura biologica in una serie di parametri matematici ( shape analysis) 
in luogo delle misure morfometriche tradizionali, consente di ridurre notevolmente la componente soggettiva 
e di sostituire o affiancare, in molti casi, il metodo tradizionale di misurazione. Nella shape analysis, la forma e 
ricondotta ad un unica computazione, che, a sua volta, riassume in se tutte le misure quali altezza, larghezza, 
curvature ed angoli. Si ottiene quindi quella che si potrebbe definire la "formula geometrica" che descrive 
I'esemplare nella sua globalita, che puo essere convenientemente trattata in analisi statistiche successive. 

II metodo e applicabile a strutture di qualunque genere, tenendo sempre presente la necessity di operare 
confronti tra punti tra loro omologhi. Alio stato attuale della tecnologia comunemente disponibile, le con- 
chiglie dei molluschi bivalvi sembrano prestarsi meglio a questo tipo di analisi, in quanto piu facilmente ri- 
conducibili ad un profilo bi-dimensionale di punti tra loro omologhi. 

In questo lavoro si confrontano i risultati della morfometria tradizionale (peso, lunghezza, altezza e spesso- 
re di ciascuna valva) e di alcuni metodi facilmente accessibili di shape analysis, basati sull'analisi di Fourier, 
applicati ad uno stesso gruppo di 39 valve destre di Goodallia triangularis. 

Sono stati utilizzati softwares differenti per trasformare le immagini digitali degli oggetti in parametri ma¬ 
tematici. Una volta acquisita I'esperienza sufficiente, il tempo necessario per passare daH'immagine digita- 
lizzata alia serie di Fourier pronta per I'elaborazione statistica e di pochi minuti per ciascuna valva. Per pas¬ 
sare dall'acquisizione delle immagini digitali al risultato finale dell'analisi statistica multivariata sono state 
impiegate complessivamente poche ore. 

II dendrogramma e la PCA eseguite sui profili delle 39 valve di G. triangularis mostrano una separazione in 
due gruppi, corrispondenti alle due forme di Goodallia triangularis: una con profilo tendenzialmente trian- 
golare e una di forma piu allungata, in accordo con precedenti conclusioni di altri autori. 

Nel caso considerato, la shape analysis e risultata notevolmente piu rapida e ha prodotto risultati piu at- 
tendibili e di piu semplice interpretazione rispetto ai metodi morfometrici tradizionali. 
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in many diverse disciplines such as palaeontology, bio¬ 
geography, ecology and palaeoecology and their appli¬ 
cations. Biological structures and shapes are often very 
complex and their morphometric description triggers at 


Introduction 

Identification of shelled species based on external mor¬ 
phology is still the basis of all subsequent elaborations 
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least four major problems: 1) the inadequacy of the use 
of a limited number of linear variables to describe a 
three-dimensional object; 2) the time-consuming job re¬ 
quired when trying to use as many linear variables as 
possible; 3) the arbitrary choice of characters to be mea¬ 
sured, that commonly leads to difficulties in comparing 
results from different researchers; 4) in small-sized, the 
ciifficulty to achieve a quantitative description of their 
morphology and the incidence of operator precision in 
the accuracy of measurements. 

Recent advances in information technology made avai¬ 
lable some user-friendly softwares for image manipula¬ 
tion and multivariate statistical analysis. In some in¬ 
stances these softwares are freely accessible from the 
web and possess a wide pool of potential users. 

The opportunity to use low-cost and friendly technolo¬ 
gies with the aim to reduce subjectivity (in morphologi¬ 
cal data acquisition) and loss of time (for sample and 
data manipulation), has been approached since the late 
eighties. It has been shown that Fourier shape analysis 
(FSA) is particularly useful for numerical description of 
diverse organisms or biological structures such as coc- 
coliths (Garratt & Swan, 1997), leaves (White et aL, 
1988), foraminifers (Healy-Williams, 1984; Belyea & 
Thunell, 1984), bivalves (Ferson et aL, 1985; Crampton, 
1996; Crampton & Maxwell, 2000), trilobites (Foote, 
1989; Cronier et aL, 1998), ostracodes (Burke et aL, 1987; 
Foster & Kaesler, 1988), blastoid echinoderms (Waters, 
1977), bryozoans (Anstey & Delmet, 1973) and fish otho- 
liths (Gauldie & Crampton, 2002; Lombarte et aL, 2006). 
FSA takes an outline contour, described as a polygon of 
digitized XY-coordinates, and decomposes it into a 
spectrum of harmonically related trigonometric curves. 
Principles of FSA have been developed into different 
transfomations. In particular, the Elliptic Fourier Analy¬ 
sis (EFA; Kuhl and Giardina, 1982) and the Fast Fourier 
Transform (FFT) have been used within DOS-based 
softwares devoted to shape analysis of fossils (Cramp¬ 
ton, 1995). 

The analysis of outlines is considered a powerful tool 
for morphometric studies, especially when homologous 
structures are lacking and consequently, landmark-ba¬ 
sed measurements do not apply (Harper, 1999). Fourier 
shape analysis converts each digitalized outline into a 
suite of numerical descriptors (the Fourier coefficients) 
that describe curves (Crampton & Maxwell, 2000). The 
coefficients are the input to multivariate statistical ana¬ 
lyses, that elaborates a numerical comparison of shapes. 
Haines & Crampton (2000) observed that EFA yelds a 
relatively large number of Fourier coefficients that are 
not computationally independent of each other. This re¬ 
dundancy would hamper and compromise the statisti¬ 
cal multivariate analysis of coefficients. To overcome 
this problem, it has been proposed another method of 
FSA called Fast Fourier Transform (FFT; Crampton, 1995). 
The aim of this paper is to suggest a friendly, easily ac¬ 
cessible and practical method for a statistical approach 
to shape analysis and to explore the sensitivity of both 
EFA and FFT by comparison with results of traditional 
132 morphometries. 


Material e methods 

Shell size is here intended as a proxy for age, even 
though the direct correlation between size and age is 
not equivalent for all specimens (Jones, 1988). 

Digital images of 39 right valves of the small bivalve 
Goodallia triangularis (Montagu, 1803) were obtained 
from a microscope Olympus BX 50 equipped with a di¬ 
gital camera Camedia C 2000z and software Camedia 
DP-soft. The digital pictures have been transformed in¬ 
to simple shell outlines (profiles), and then into a series 
of mathematical parameters to be analysed by multiva¬ 
riate statistics (Fig. 1). 

Presently, the entire procedure requires four different 
soft-wares. A common graphic software to extract the 
profile; a software for image manipulation and photo 
touch that can produce the co-ordinates of each pixel of 
the profile; a common work-sheet to store the co-ordi¬ 
nates and their Fourier transform; "PAST", a free soft¬ 
ware for statistical analyses in biology and palaeonto¬ 
logy (Hammer et aL, 2005) allowing to standardize data 
and to perform multivariate statistical analysis of profi¬ 
les (similarity matrix, dendrogram and Principal Com¬ 
ponent Analysis). FFT coefficients have been obtained 
with the programs Hmatch and Hangle (Crampton & 
Haines, 1996) and subsequently elaborated with PAST 
in order to perform the analysis of similarity and the 
PCA ordination. 

The valves have been randomly extracted from the 250 
right valves on which a morphometric analysis based 
on weight (as a proxy of shell thickness), length, height 
and width (as synthetic descriptors of shape) of each 
valve was already performed (Basso, 1992). This ap¬ 
proach allows the comparison of the results obtained 
with the traditional method and with the digital-profi¬ 
les method. Procedural details are listed in appendix. 

Results 

The multivariate statistical analysis of outlines of 39 G. 
triangularis right valves produced graphic results (den¬ 
drogram and PCA scatter plots, Figs. 2-5) in which each 
point represents one right valve profile. The statistical 
analysis produced different results for EFA- and FFT-ba- 
sed shape analysis. Some representative outlines are re¬ 
ported over the corresponding number in the PCA plot. 

Elliptic Fourier Shape Analysis 

The dendrogram (Euclidean distance, paired groups. 
Fig. 2) and the PCA (Fig. 3) based on the 39 profiles of G. 
triangularis right valves show three clusters (A to C). 
Cluster A is made of ventrally rounded valves (G. trian¬ 
gularis forma subtrigona), cluster B contains valves with 
more triangular shape (G. triangularis forma typica) (Fig. 
2). Cluster C contains valves with intermediate shapes 
and the shell outlines at its border could be included in 
the first and second clusters. In the third cluster, outlines 
no. 37 and 5 correspond to intermediate shapes. Results 
of PCA (Fig. 3) show that the first two components re¬ 
presents about 91% of the total variability. 


Fig. 1. Schematic procedure to perform 
multivariate statistical analysis of shell out¬ 
lines. Details in the Appendix. 

Fig. 1. Schema della procedura per esegui- 
re un'analisi statistica multivariata su profili 
di conchiglie. Dettagli in Appendice. 



Fast Fourier Transform 

The dendrogram based on coefficients produced by FFT 
(Euclidean similarity, paired groups, Fig. 4) shows only 2 
clusters (A and B), represented also in the PC A scatter plot 
(Fig. 5). Cluster A includes G. triangularis form subtrigona 
and cluster B encompasses specimens of G. triangularis 
form typica. The shells with the maximum lenght/height 


ratio are concentrated between values of 0.02 and 0.01 of 
the first component and between 0 and 0.04 of the second 
component. Tire intermediate shell outlines that were 
plotted in cluster C (EFA plot) are here splitted into clu¬ 
sters A and B. In particular, outlines no. 5, 12 and 37 are 
assigned to cluster B, close to the boundary with cluster 
A. As a result, FFT produces a sharp separation of the out¬ 
lines into clusters. 


Fig. 2. Dendrogram produced by hierar¬ 
chical clustering of trigonometric coeffi¬ 
cients (corresponding to the outlines) ob¬ 
tained with EFA. Clusters A to C are obtai¬ 
ned at the - 0,015 level of Euclidean dis¬ 
tance. Clustering calculated on unweigh¬ 
ted pair-group average (UPGMA). 

Fig. 2. Dendrogramma prodotto dal clu¬ 
stering gerarchico agglomerativo dei coef¬ 
ficient trigonometrici (corrispondenti ai 
profili) ottenuti con I'Analisi Ellittica di Fou¬ 
rier (EFA). I gruppi da A a C si separano al¬ 
ia distanza euclidea = -0,015. II clustering 
e calcolato sul legame medio (UPGMA). 


LOcsjT-^cooih-cDooooouO'srr^co^— cnj ^ Or-o r^ocMtncoco^cooo 

^ CM CO ^ CO ^ T- CM ^ CM CM CM CsJ CN CM CN CO CO -sf- CO CO CM OO CD CO CO T-LO CO 'I-^ CO ^ CM CO CO 



133 


Canonical morphometry versus statistical treatment of outlines through Fourier shape analysis: an empirical comparison 










































































































Francesco Brusoni & Daniela Basso 


134 



Discussion 

In order to study the morphometric variability of Goo- 
dallia triangularis occurring in a Mediterranean shell as¬ 
semblage (= thanatocoenosis). Basso (1992) faced the dif¬ 
ficulty to measure and deal with valves that rarely ex¬ 
ceeded 3.3 mm in length. It was required a time-consu¬ 
ming job and the use of high-precision calibres and ba¬ 
lance. Subsequently bivariate and multivariate statisti¬ 
cal analyses of the measured variables allowed the reco- 
gnitiion of two ecomorphotypes of G. triaugularis, one 
more triangular, representing the Atlantic, cold-water 
form and now occurring in the relict sediment of the 
Mediterranean shelf-break; the other more elongate, 
which is the Mediterranean living form (form snbtri- 
gona) (Fig. 6). These conclusions were later confirmed by 
Giribet & Penas (1999) in the framework of the revision 
of the genus Goodallia. The occurrence of intermediate 
specimens, as observed by Basso (1992), prevents the se¬ 
paration of these eco-morphotypes into different spe¬ 
cies. In the same paper (Basso, 1992), most of these inter¬ 
mediate specimens, interpreted as juveniles on the basis 
of their very small size, were positioned at the left end 
of the PC A scatter plot (Basso, 1992, Fig. 6). The largest 
specimens were spread along the right end of the PCA 
plot, with the longest and relatively light valves at the 
bottom right (form subtrigoua), and the highest, relati¬ 
vely heavy valves at the top right of the plot (form 
typica; Basso, 1992, Fig. 6). 

The complete outline of the shells of the two subspecies 
of Goodallia triaugularis summarizes the general shape, 
the slope of the margins and the roundness, and the 
more or less equilaterality of the bivalve on the basis of 
the position of the umbo (orthocline in the subspecies 


Fig. 3. PCA scatter plot based on the same 
coefficients (EFA) and Euclidean distance of 
Fig. 2. The A to C clusters are indicated. 

Fig. 3. Grafico PCA basato sugli stessi 
coefficient! (EFA) e distanze euclidee di Fig. 
2. Sono indicati i gruppi da A a C. 

typica, prosocline in the subtrigoua). Unlike "traditional" 
morphometries, shape analysis considers only profiles 
and disregards dimensions. Therefore, unlike Basso 
(1992) our results show juvenile shells within all three 
clusters of EFA-based dendrogram, to indicate a possi¬ 
ble early divergence in shell morphology. 

When a convenient number of specimens is available, it 
is possible to study the reciprocal position of the 
"cloud" of points in the PCA plot, along possible gra¬ 
dients of shape variability. 

It is of a paramount importance to verify that the samples 
are photographed with exactly the same position. Al¬ 
though PAST contains a standardization procedure allo¬ 
wing to deal with some variability in sample orientation, 
it is better to homogenize the orientation in order to avoid 
confusion or subjective bias in the subsequent phases. 
Moreover, it must be considered that possible displace¬ 
ments of a bivalve shell are in one plane only, and that is 
quite easily managed with PAST standardization. Major 
problems are expected to arise from gastropod shell orien¬ 
tation, since rotation along the long axis of the shell or til¬ 
ting will result in a different profile on the plane, and thus 
the same shell variably oriented will be recognized as a 
different object. This kind of biased image acquisition can¬ 
not be solved by PAST standardization (Hammer et al, 
2005). 

Most difficulties arise in establishing an homology bet¬ 
ween points of profiles; points must be philogenetically 
homologous in order to perform a meaningful statisti¬ 
cal analysis. (Bookstein et al, 1982; Read & Lestrel, 1986; 
Temple, 1992). In many cases it is possible to identify 
homologous points in a biological specimen. These 
points are taken as landmarks in traditional morpho¬ 
metry Thousands of literature examples cover very di- 





Fig. 4. Dendrogram produced by hierar¬ 
chical clustering of trigonometric coeffi¬ 
cients (corresponding to the outlines) ob¬ 
tained with FFT. Clusters A and B are ob¬ 
tained at the - 0,017 level of Euclidean 
distance. Clustering calculated on un¬ 
weighted pair-group average (UPGMA). 

Fig. 4. Dendrogramma prodotto dal clu¬ 
stering gerarchico agglomerativo dei coeffi¬ 
cient! trigonometric'! (corrispondenti ai pro- 
fili) ottenuti con Trasformata Veloce di 
Fourier (FFT). I gruppi A e B si separano alia 
distanza euclidea di - 0,017. II clustering e 
calcolato sul legame medio (UPGMA). 



verse subjects such as trilobites, molluscan shells or the 
biometry of human skull (for example: Hughes & 
Chapman, 1995; Lacroce & Repetto, 2004). In other in¬ 
stances these points that are easily identified from a 
biological point of view, are difficult to locate with pre¬ 
cision in the practice. This can be the case of the bivalve 
umbo that corresponds to a small area on the shell ra¬ 


ther than a precise point (Crampton, 1995; Crampton & 
Haines, 1996). This problem represent a possible source 
of uncertainty in manual measurements, and is negligi¬ 
ble only for large-sized specimens. Similarly, it is diffi¬ 
cult to locate homologous points on the spire of gastro¬ 
pods, with the exception of the apex (protoconch). The 
use of a complete outline, as in Fourier shape analysis, 


Fig. 5. PC A scatter plot based on the same 
coefficients (FFT) and Euclidean distance of 
Fig. 4. The A and B clusters are indicated 
(from Basso, 1992). 

Fig. 5. Grafico PCA basato sugli stessi 
coefficient! (FFT) e distanze euclidee della 
Fig. 4, Sono indicati i gruppi A e B (da Bas¬ 
so, 1992). 
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corresponds to the use of all homologous points at the 
same time, therefore the identification of selected land¬ 
marks is no more necessary (Foote, 1989). However, all 
specimen outlines must be obtained with the same spe¬ 
cimen orientation, in order to allow comparison among 
series of homologous points. 

Haines & Crampton (2000) observed that EFA yields a 
relatively large number of Fourier coefficients that are 
not computationally independent of each other. This re¬ 
dundancy would hamper and compromise the statisti¬ 
cal multivariate analysis of coefficients. Actually, FFT of 
our Goodallia outlines produces the sharpest separation 
of the shells into clusters, which is an obvious advantage 
for interpretation. EFA results are comparably less clear. 
Interestingly, the same outlines (15 and 36) are at the op¬ 
posite ends along component 1 in both PCA, while outli¬ 
nes 23 and 3, that are opposite along component 2 in EFA 
(Fig. 3), are also split by component 1 in FFT (Fig. 5). 

The procedure used here is much faster and simple in 
comparison with traditional morphometries and allows 
an easier graphic interpretation of results. Despite the 
use of several different softwares, the entire procedure 
from digital photographs to dendrogram and PCA plots 


Fig. 6. A. Principal component analysis 
for L = length; H = height; W = width; Wt 
= weight. B. component scores plot (from 
Basso, 1992). 

Fig. 6. A. Analisi dei component! - princi¬ 
pal (PCA) per L = lunghezza; H = altezza; 
W = spessore; Wt = peso. B. Grafico dei 
component scores (da Basso, 1992). 

required few hours (Fig. 1). With some practice, only 
few minutes are necessary to produce the Fourier series 
from the digital picture. 

Conclusions 

This paper is aimed to compare the performance of tradi¬ 
tional morphometries and of FSA in the description of 
small bivalves. In agreement with Davoli (1982), the best 
practice of a morphologic-comparative study should be 
the use of observations and measurements of all morpho¬ 
logical variables, avoiding the use of few, arbitrarily cho¬ 
sen characters of classification. In the case of Goodallia , the 
very small size and the fragility of shells prevent the col¬ 
lection of further linear measurements in addition to shell 
length, height and width considered by Basso (1992). 
Moreover, some morphologic features are traditionally 
defined in a qualitative way (such as the symmetry of the 
bivalve shell and its roundness), preventing their quanti¬ 
fication, objective definition and statistical treatment. 

The first advantage of FSA is the absence of assump¬ 
tions on characters to be selected and measured. The 
transformation of the outline of a biological structure 








into a suite of mathematical parameters reduces subjec¬ 
tivity since outlines are the synthesis of many measure¬ 
ments (height, length, width, slopes, convexity, etc.). It 
can be argued that some features cannot be detected by 
FSA of outer side of bivalve shells, such as internal 
structures, shell thickness or ornamentations. However, 
most internal features of bivalves correspond to major 
supraspecific diagnostic characters, such as the type of 
pallial line (integer or with a sinus) and the hinge struc¬ 
ture, that are efficiently and quite rapidly described in 
the traditional qualitative approach. Moreover, thick or¬ 
namentation such coarse ridges and costae, spines etc. 
are expected to be detectable along the outline. 
Multivariate statistical elaboration had analogous bearing 
in respect to the traditional morphometry, leading to an 
easy separation of the two forms of Goodallia with a much 
faster and less subjective procedure. The same results, ob¬ 
tained with the canonical morphometry based on linear 
measurements (Basso, 1992), required electronic instru¬ 
ments for precision mechanics and a time-consuming 
specimen manipulation. The second advantage of FSA is 
the objectiveness and the rapidity of data acquisition. 

In the comparison of the PCA results based on the two 
different Fourier transformations, FFT produces the 
sharpest separation of the outlines into clusters while 
EFA results are comparably less clear (Figs 2- 3). 

At the present state of the art, FSA is not expected to 
substitute all canonical morphometric methods. A criti¬ 
cal point remains the choice of topics that can be succes¬ 
sfully afforded through the comparison of bi-dimensio- 
nal profiles of homologous points, as most bivalves are. 
Gastropods could be eligible for FSA, provided that a 
constant orientation of the specimens is ensured. 
Among bivalves, the analysis of similarity among diffe¬ 
rent species of the same genus or among infraspecific 
taxa appears particularly promising. 

The method described in this paper involves the use of 
different softwares. Work is in progress to merge and in¬ 
tegrate them into one single software that would be a step 
forwards the divulgation of such a powerful technique. 

Appendix 

Procedure : 

1) Digital photographs 

Specimens must be complete and broken margins must 
be avoided, since they will be interpreted by FSA as a 
dramatic change of outline. The color of the specimen 
background should be selected to contrast the shell outli¬ 
ne and enhance separation while working at pixel level. 

2) Obtaining outlines 

All digital pictures to be dealt with must have standard 
size (pixels) in order to obtain black outlines, composed 
of the same number of points, on a white background 
(colours can be inverted). This condition must be sati¬ 
sfied when entering the PAST procedures for similarity 
analysis. Common graphic softwares can be used to ob¬ 
tain the simple shell outlines and "clean" them from 
possible "noise" pixel generated by dust, shadows etc. 


3) Transforming the outline graphic file into a suite 
of co-ordinates (X, Y) 

Most common softwares will "read" the outline the occi¬ 
dental way: from left to right, from top to bottom, as a 
written text. On the contrary, the outline transformation 
into co-ordinates requires to read it counterclockwise, each 
pixel only once. Therefore it is necessary to divide each 
valve profile into half (right and left) the will be dealt with 
separately. Successively the two halves will be converted 
in ASCII (.dat) file containing the points co-ordinates. 

The ASCII file contains the description of the profile given 
by the co-ordinates of all the composing points. Opening 
the electronic sheet of each half profile we will find three 
columns: the first for abscissas, the second for ordinates, 
the third for the RGB colour of each point. Using filters it 
is possible to separate the profile points from the back¬ 
ground simply selecting the colour (black points = RGB 0). 
After this stage a new file should be created containing 
only the two columns of co-ordinates. 

The same procedure is repeated for the other half profile 
that must be put upside down, in order to allow joining 
the first half reading the whole profile counter clockwi¬ 
se. This way the table will show the co-ordinates of 
points of the profile starting from the uppermost point 
of the left half-profile, down to the lowermost point of 
the same left half-profile, up in a counter clockwise sen¬ 
se from the lowermost point of the right half-profile to 
the uppermost point of the same right half-profile. 

All couples of co-ordinates (the first for abscissas, the 
second for ordinates) pertaining to all outlines should 
be saved in the same file. 

Example: 

XI a Y1 a XlbYlb . where each letter a, b, ... 

X2 a Y2 a X2 b Y2 b identifies an outline 

X3 a Y3 a X3 b Y3 b 


4) Transformation of co-ordinates into coefficients 
Two different pathways are followed depending on the 
selected type of Fourier transform: Elliptic Fourier 
Analysis or Fast Fourier Transform. 

Elliptic Fourier Analysis (EFA) 

PAST programs 

The co-ordinates file is copied into PAST (Hammer et al., 

2005). It is now necessary to transpose the matrix of co¬ 
ordinates on a line, alternating abscissas and ordinates: 

Example: 

Xla Yla X2 a Y2aX3aY3a . 

X1Z? Ylb X2b Y2bX3bX3b . 

This can be obtained with the PAST function " GROUPED 
COLUMNS TO MULTIVAR", (edit menu) indicating two di¬ 
mensions. Then the transpose function can apply. 

In order to standardize profiles in size and rotation on the 
plane, it is necessary to use the PAST function PROCRUSTES. 

This function can be lengthy but is necessary to give pro¬ 
files the same orientation. Now profiles are converted 137 
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with the function Elliptic Fourier shape analysis in a se¬ 
quence of trigonometric coefficients (Ferson et ah, 1985). 

Fast Fourier Transform (FFT) 

HANGLE, HMATCH programs 

The co-ordinates file is the input for the program Han- 
GLE (Crampton & Haines, 1996) that output the Fourier 
coefficients. Successively, the Hmatch program normali¬ 
zes the Fourier coefficients of each outline for orienta¬ 
tion and starting position of the trace. 

5) Multivariate statistical analysis 
PAST can also be used for classification and ordination 
of coefficients obtained with EFA or FFT (cluster analy¬ 
sis and Principal Components Analysis = PC A). In this 
paper we used Euclidean distance and unweighted 
pair-group average (UPGMA). PAST has other func¬ 
tions devoted to improve the interpretation of PCA or¬ 
dination, such as the Minimum spanning tree or the 3D 
option for the multidimensional disposition and rela¬ 
tionship of points in the plot. 
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